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Abstract: 

We compute non-perturbatively the average up-down and strange quark masses from the 
large momentum (short-distance) behaviour of the quark propagator in the Landau gauge. 
This method, which has never been apphed so far, does not require the exphcit calcu- 
lation of the quark mass renormalization constant. Calculations were performed in the 
quenched approximation, by using (9(a)-improved Wilson fermions. The main results of 
this study are (2GeV) = 5.8(6) MeV and mf{2GeV) = 136(11) MeV. Using the 
relations between different schemes, obtained from the available four-loop anomalous di- 
mensions, we also find mf^^ = 7.6(8) MeV and mf'^^ = 177(14) MeV, and the MS masses, 
mp(2GeV) = 4.8(5) MeV and m^{2GeV) = 111(9) MeV. 



PACS numbers: 12.15Ff, 11. 15. Ha, 12.38.Gc. 



1 Introduction 



Determination of quark masses is becoming one of the most intensive field of investigation 
in lattice QCD The accuracy of the predictions is significantly improving mainly 

because of two recent theoretical developments: 



• Non-perturbative renormalization procedures have been introduced ||15| |I6| in order 
to remove the systematic uncertainties coming from the truncation of perturbative 
series in calculation of the relevant renormalization constants. These procedures 
also provide an appropriate non-perturbative, short distance definition of the quark 
masses either in the so-called RI-MOM or Schrodinger functional schemes. The 
relation between the mass in the RI-MOM scheme (which will be used in this study) 
and in the MS scheme [], or the renormalization group invariant mass, is known at 
next-to-next-to-leading order (N^LO) [l^, and very recently even at N^LO [|l^, in 
continuum perturbation theory. 

• A second important theoretical progress is the reduction of finite cut-off (C(a)) effects 
obtained by improving the lattice fermion action and operators. The perturbative 
procedure, proposed in Refs. [|19|, has been recently extended to a fully non- 
perturbative (9(a)-improvement by ALPHA collaboration [^], so that the remaining 
discretization errors are only of 0{a'^). 

In the past year, several independent lattice determinations of the light quark masses 0- 
iTI have been presented, adopting both non-perturbative renormalization procedures and 



non-perturbative improvement. 

The two standard definitions of the lattice quark masses are based on the vector (VWI) 



and the axial- vector (AWI) chiral Ward identities ||22|. The VWI relates the bare quark 
mass to the value of the Wilson hopping parameter, 2am = — 1/Kcrit)- With this 
definition, one can easily show that the mass renormalization constant is ZmifJ^) = Zg (/i). 
The definition based on the AWI is 2am = {a\d^A^\(3) / {a\ P\P), where d^A^ and P are the 
divergence of the (improved) axial vector current and the pseudoscalar density, respectively. 
In this case, Zf^^n) = ZA/Zp{fi). 

In this paper, in order to calculate the renormalized quark mass, we adopt a new method 
based on the study of the large-p^ behaviour of the renormalized quark propagator. The 
method is based on the idea that at large Euclidean momenta it is possible to match 
lattice and continuum correlators by requiring the vanishing of chirality violating form 
factors This procedure is justified by the following two observations. The first 

is that at large momenta the renormalized perturbation theory becomes chirally invariant 
(explicit chiral symmetry breaking effects induced by the regularization are reabsorbed 
by imposing the validity of the chiral Ward identities, while violations from the non- 
vanishing quark masses, disappear at large momenta). The second observation is that the 

^ The conversion of the resuhs to the MS scheme is only necessary for comparison with other calculations 
for which the quark masses are renormalized perturbatively. Otherwise, the method of Ref. [ p^ allows, 
in principle, to obtain the renormalized quark masses for the RI-MOM scheme in a completely non- 
perturbative way. 
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contributions due to the spontaneous breaking of chiral symmetry, which are absent in 
perturbation theory, die off at large momenta. Thus, both effects decrease as we go deeper 
into the Euchdean region. 

The simplest application of this idea is the possibility of relating the quark mass to the 
renormalized quark propagator 



S[P) = + T" 



since, at large p^, we expect p6 



a2{p') ~ m + {qq)^ + 0{l/p^) . 



(2) 



In particular, the quark mass renormalized in the RI-MOM scheme, can be directly ex- 
tracted from the quark propagator renormalized in the same scheme by using 



m^^H/i) = -Tr S-\p■^i) 



(3) 



where the trace is over both color and spin indices 0. S{j>] n) is the (improved) quark 
propagator renormalized at some scale /i. Since the quark propagator is a gauge dependent 
quantity, the definition of the RI-MOM mass also depends on the gauge. 

At large momenta and up to discretization errors, Eq. (|]) is equivalent to the definition 
of the quark mass based on the AWL Chiral symmetry provides a relation between the 
inverse quark propagator, S{p] and the amputated Green function of the pseudoscalar 
density, A5(p; /i), computed between external (off-shell) quark states of equal momenta p. 
The AWI then reads: 



2mg(/i) A5(p;/i) = 75S ^(p;/i) + 5 ^(p;/i)75 



(4) 



All quantities in Eq. (|) are assumed to be renormalized (and improved) in the same 
scheme, and at the same scale /z. In the RI-MOM scheme (and in a fixed gauge), the 
Green function h^{p] fi) satisfies the following renormalization condition [p!5[] : 



In. 

12 



ib^bip; /i) 



(5) 



By tracing both sides of Eq. (H) with 75 and by using Eq. (|[), the relation (|) is readily de- 
rived. Note that if /i in Eq. @ is not chosen in the perturbative region, i.e. fi ^ ^qcd, the 
definition of the quark mass will be affected by non-perturbative, chirally-breaking contri- 
butions proportional to the quark condensate and higher-dimensional operators, appearing 
in higher power corrections (oc 1/p^"). 

The advantage in determining the masses from the quark propagators is that it is not 
necessary to calculate explicitly the mass renormalization constants {i.e. Zs{fi) or Zp{fi)). 

^We note in passing, that by using this method we were not able to extract the value of the quark 
condensate (||). 
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This is merely a consequence of the fact that the renormahzed quark propagator is directly 
expressed in terms of the renormalized quark mass. Unlike in the case of the VWI, the 
critical value of the hopping parameter, Hcrit, is also not needed (which is the advantage 
inherent to the use of the AWI). There is, however, one renormalization constant for any 
quark- mass definition: in our case, using Eq. (^), this is the quark- field renormalization 
constant Zg. In the RI-MOM scheme Zg is fixed by the following renormalization condition: 



^Tr 

48 



48 



7/.- 



d S-\p) 



(6) 



In summary, Eqs. @ and (§) are all we need to extract quark masses from propagators. 
The procedure becomes rather complicated, however, if we want to extend it to the non- 
perturbatively improved case. The drawback with the method of Ref. |jl5| is that the 
improvement program (which was initially carried out for on-shell quantities) must be 
extended to off-shell Green functions on non-gauge invariant states and involves additional 
counterterms for a full 0{a) improvement The strategy followed in this case will be 
illustrated in detail in Sec. 

In this paper, Eqs. (H) and @ have been used to compute the average up-down and the 
strange quark masses, by performing a lattice QCD calculation in the quenched approxima- 
tion. We use the non-perturbatively improved action ||2T|, and improve the quark propaga- 
tor in the chiral limit. The (9(a)-improvement procedure for off-shell (gauge non-invariant) 
quantities has been discussed in Ref. p3[. Since that paper has not been pubhshed yet. 



we will describe here in some detail the specific case of the quark propagator. For techni- 
cal reasons, which are related to the mixing with non gauge-invariant higher- dimensional 
operators (see below), we are not able to improve the propagator out of the chiral limit 0. 
Therefore, our determination of the quark masses is affected by 0{gQam) systematic errors. 
Since the value of the inverse lattice spacing in this simulation is ~ 2.72 GeV, these 
errors are expected to be negligible for the strange and the light quark masses. 

We conclude this section by summarizing the main results of this paper. From the 
study of the quark propagator, we extract the (quenched) light and strange quark masses 
in the RI-MOM scheme: 



mf\2 GeV) = 5.8(6) MeV , (2 GeV) = 136(11) MeV 



RI/ 



(7) 



These results are in very good agreement with those of Ref. 0, namely (2GeV) = 
138(15) MeV, and mf\2GeV) = 5.6(5) MeV. 

Using the N^LO perturbative formulae of Ref. [0, we obtain the renormalization group 
invariant quark masses 



m 



RGI 



7.6(8) MeV 



m 



RGI 



177(14) MeV 



where the renormalization group invariant quark mass is defined according to the conven- 

i.e. 



tion usually adopted in perturbative calculations |]T^ Il8|, BS 



m 



RGI 



hm m,(/i) (a,(/i))-^™ 



/I— >oo 



(9) 



^ The general problem of the improvement out of the chiral limit has not been solved yet, although 
several interesting proposals exist |23, kJ, 27, 2q|. 
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with 7^"* and /3o, being scheme independent Q 

Finally, the quark masses in the MS-scheme read 

NLO N^LO N^'LO 

mp(2GeV) = {5.2(5); 4.9(5); 4.8(5)} MeV , 

mf^(2GeV) = {120(9); 114(9); 111(9)} MeV , (11) 

where the numbers within the curly brackets are obtained after converting the RI-MOM 
results to the MS one to NLO, N^LO and N^LO accuracy, respectively. The details on 
anomalous dimensions and beta function are listed in the appendix. 

We note that in most of the phenomenological applications, for example with QCD 
sum rules, the theoretical expressions are only known to the NLO and, for consistency, the 
quark masses at the same accuracy should be used. 

Similarly, we stress that lattice calculations of quark masses, in which the mass renor- 
malization constants has been determined by using (one-loop) perturbation theory, should 
be compared with our NLO results of Eq. (|ll|), since they have been derived at the same 
order of accuracy. 

Preliminary results obtained with the method discussed in this paper were already 



presented at the "Lattice 99" Conference |30 . 



2 Improved quark propagator 

The general problem of improving gauge non-invariant, off-shell, correlation functions has 
been studied in Ref. [^. Since this paper is still unpublished, in this section we discuss 
in some detail the non-perturbative improvement of the lattice quark propagator to 0{a). 

2.1 The subtracted quark propagator 



The original idea of improvement (later developed in Ref. for gauge theories) 
consists in adding, to both the action and operators, a complete set of higher- dimensional 
( ^'irrelevant") operators, the coefficients of which are tuned as to cancel finite cut-off effects 
(to a desired order of lattice spacing). Specifically, the improvement of the Wilson action 



to 0{a) is achieved by adding a set of dimension- five operators ^1 



n „ 

S = Sw + aJ2ci d^x Of='^ (x) , (12) 



''Note that in Ref. ||T^, another convention has been used: 



'2^0 



m^Gi = Urn TO,(^) ( ^as{^Ji) ] . (10) 
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allowed by gauge invariance and discrete lattice symmetries, namely 

i 1 

03 = m^qq, = m q (p + mo)q , O5 = q(p+mofq, (13) 



where m is the bare subtracted mass and p + mo the bare Dirac operator appearing in Sw- 
The operators O2 and O3 can be reabsorbed in the definition of the bare strong coupling 
and the quark mass. A major simplification comes from the restriction of improvement 
to physical amplitudes, for which the equations of motion can be used. In this way, one 
is left with one (Clover) counterterm only, Oi, the coefficient of which was computed 
non-perturbatively in Ref. pl| . 

The equations of motion cannot be used to improve off-shell quantities (such as the 
quark propagator): in this case, one must also consider the operators 04^5. As for operators 
which are BRST allowed but not gauge invariant, only the BRST variation of c°'d^A^ 
(where and are the anti-ghost and gauge fields, respectively) may contribute. This 



term, however, can be absorbed into a redefinition of the gauge-fixing parameters |23 



Besides the terms in the action considered above, in Ref. it has been also shown that 



for the quark propagator, there is another operator, which is not BRST invariant but 
may contribute to off-shell correlation functions (because its presence is not excluded by 
Slavnov- Taylor identities). The effects of 04 5 and of this extra operator can be eliminated 
with a simple redefinition of the quark field: 

q{x) =Zf^ -^'^{l+bqma){l + ac'^^p + m^) + acNci^} q{x) . (14) 

We now discuss how the unknown coefficients Zq {Zq^^'^ = zf^ ^^"^{1 + h^ma)), c'^ 
(corresponding to the coefficient of the operator O5) and cngi, present in Eq. (|I^), can be 
determined from the analysis of the lattice bare quark propagator, Sl{p). 

From Eq. (p^lf ), it follows that the relation between Sl{p) and the improved, renormal- 
ized quark propagator, S{p), constructed in terms of the quark fields, q and g respectively, 
has the form: 

SLip) = il-2ac^Gii^)ZgSip)-2ac'q . (15) 

In this equation, it is convenient to express the renormalized quark propagator, S{p), in 
terms of the two invariant scalar form factors, 0^1 (p^) and 5^2(p^); defined in Eq. (|I|). For 
further use, we remark that at large p"^, up to power-suppressed (~ 1/p^) and logarithmic 
corrections, ?i(p^) — 1 and cr2(p^) — m, where m is the renormalized quark mass. After 
substituting (P into ([l5|) , one finds: 

(TiLip^) = ^ Tr [-2|^^l(p)] = Zg (a,{p^) - 2a cngi Mp^)) (16) 



— Tr 



p2 12 



Sl{p) = -2a Cq + 2ac^GiZqdi{p ) + Zq (17) 
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where crii^2L(p^) are the analog of ai_2(p^) for the lattice bare propagator. 

Using Eqs. (|T6|) and (|I7D , the coefficients and ^^cngi can be determined as follows: 
at large and in the chiral limit, since ^2 ~ — > 0, it is in principle possible to 
separate and ZqC-^Qi using the dependence of ai; the overall renormalization constant 
Zg, including its 0{a) mass dependence, can then be determined by combining Eq. (|16]) 
with the renormalization condition (^. In practice, however, this procedure is very difficult 
to implement. The reason is that the logarithmic p^-dependence of ai(p^), entering the 
right-hand side of Eq. (jlT]), is very mild. This is especially true in the Landau gauge, where 
it starts at order a1 in perturbation theory. Thus, it is very hard (if not impossible) to 
disentangle the contributions coming from the two coefficients Q c'^ and cngi- 

Let us return to the large-p^ behaviour of the lattice quark propagator, in the limit 
where power and logarithmic corrections can be neglected. In this limit, Eqs. (|TBp and 
(p!7|) become: 

o"il(p^) ^ Zq{l- 2a Cngi "^j = Zq 



and 



2ac'+2ac^GiZq = -2aZ' (19) 



In the large-p^ region, by using Eq. (|1^), it is then possible to determine an "effective" 
renormalization constant, Zq, which reduces to Zq in the chiral limit. Moreover, in terms 
of the coefficient c^, computed through Eq. (|19|), we can define a "subtracted" quark 
propagator, S{p), as: 

S{p) = Z-^{SL{p) + 2ac'^. (20) 

If the coefficient cngi were equal to zero, the subtracted propagator S{p) would correspond 
to the improved, renormalized quark propagator, S{p), as can be seen from Eq. ([15|) . In 
the presence of cngi, however, S{p) and S{p) differ by terms of 0(cngi am), up to (small) 
logarithmic corrections. We conclude that, by following the procedure outlined above, we 
are able to exactly improve the quark propagator in the chiral limit. Out of the chiral 
limit, since cngi is of O^Qq) in perturbation theory, the propagator is affected by O {gharri) 
discretization errors. In the range of quark masses considered in this paper, these terms 
are expected to be smaller than other statistical and systematic uncertainties. They may 
be important, however, in the calculation of heavy quark masses. 



2.2 Practical implementation 

Let us now discuss how the improvement procedure for the lattice quark propagator works 
in practice. As a preliminary (and instructive) step, we consider the inverse unsubtracted 
lattice propagator expressed in terms of the usual form factors, Si^ and T,2l'- 

Sl\p) = -if^idp') + ^2l{v') . (21) 

promising way to compute separately and cngi non-perturbatively is from the study of the 
quark-gluon vertex. 
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T^iL is special in that its p^-behaviour is protected by the VWI, T^nijP') ~ const. = Zy, as 
can be seen in Fig. |l] (empty circles) 0. By using Eqs. (|I|) and (0), in the large-p^ limit 
one finds 



2 

2 



1 + p 

\p^criL[p- 



-1 



Sr2(p^) [ 1 + 0(aV) ] . (22) 



The effect of the 0{a?) term in Eq. (^) is important at large p^, as shown in Fig. |T|: 
Sil(p^) (empty circles) is indeed flat, whereas <J^l{p'^) (empty squares) is not. The reason 
is the presence of the contact term in the factor 

1 + - 1+ - 1 + , (23) 

at large p^. For c'^ = 0, instead, the factor 1 + z'^p'^ ~ 1 + m? /p^ 1. Thus the dangerous 
0{a^) terms of Eq. ([2^ ) jeopardizes the p^ behaviour that di should have in the continuum. 

The following procedure has been adopted to get rid of the 0{a?) corrections induced by 
the contact terms []. Motivated by Eq. (0), we flrst multiply the original lattice propagator 
by the overall factor, (1 + z^p^). Since the action is only (9(a)-improved, this subtraction 
is formally irrelevant. The multiplication, however, removes the dependence of (yiL{p'^) on 
p^, coming from the O(a^) effect discussed above. This can be seen in Fig. 0by comparing 
the squares and fllled circles. It is important to add that z (for each k) has been flxed by 
fltting to a constant the ratio 

in the large p^-region. The data and the flt are displayed in Fig. |1] for k = 0.1344. The 
numerical values for z will be given in the next section. The multiplication by (1 + z^p^) 
flattens the dependence of both and a2L- This effect is particularly pronounced for 
(JiL, as can be seen in Fig. |l]. 

We now describe the procedure followed to remove the constant contact term c^, in 
Eq. (^). After the multiplication of the propagator by (1 + z^p'^), we flt c'"2l(p^) to the 
form expected from the OPE, namely (up to logarithmic corrections) 

f!HiM=^+^ + CaV, (25) 

in the region of large Euclidean momenta where the OPE applies. We have used the interval 
0.5 < a^p^ < 2.0, which at /? = 6.2 corresponds to 2 GeV^ < < 15 GeV^ in physical units. 
By comparison with Eq. (|19]), the coefficient A is evidently A = —2aCg. From the fit of our 

^Throughout this paper, we adopt a continuum notation in which ap^ stands for sin(ap^). Thus, for 
instance, a^p^ corresponds to sin^(ap^) and a;^ is equal to 7^ sin(ap^). 

^ Although not mentioned before, it is clear that there are other 0{a^) effects, besides those induced by 
the 0{a) contact term, resulting in the factor (1 + z^p^). These, however, are found to be much smaller, 
see below. 
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Figure 1: The left Ggure shows the effect of the correction due to the factor (1 + z^p^). 
Empty symbols denote the quantities extracted directly from the lattice quark propagator. 
The filled circles denote the effect of the correction. Note that the error bars are smaller 
than the used symbols. The value of the constant z as obtained from a fit to our data in 
the large -region is shown in the right figure. Both figures correspond to n = 0.1344. 



data to Eq. (|25D , we obtain A = 0.617(11), to be compared to A^^^^^ = 0.573, as computed 
in one-loop (boosted) perturbation theory ^ According to Eq. (^, the parameter B is 
proportional to the quark mass, so that it is expected to vanish in the chiral limit. The value 
that we obtain, B = 0.003(6), is well consistent with expectations. It clearly demonstrates 
that the contribution of the term proportional to the quark condensate is negligible in the 
range of momenta chosen for the fit. This point was recently questioned in Ref. From 
our fit this contribution appears to be completely negligible for > 3GeV^, as expected 
from OPE, if (qq) ~ ^qcd- Finally, we also obtain C = 0.022(4), for the parameter which 
contains the information on residual O(a^) effects. Note that without the correcting factor 
(1 + z'^p'^) we would have obtained C = —0.070(3), larger than the above results. This 
shows that the residual O(a^) effects are smaller than those induced by the contact terms. 
To summarize, we subtract the unwanted discretization effects using 

S{p) = 1 5i(/) + ^ = Sl{p) (l + zY) -a- CaV . (26) 

The resulting improved propagator, S{p), exhibits a good chiral behaviour at large p"^ 
and its inverse has small O(a^) corrections. S2(j5^) is expected to be a slowly varying 
function of the momentum at large p^, with a logarithmic p^-dependence governed by the 
quark mass anomalous dimension. As shown in Fig. 0, the bare form factor (S2l(p^)) 
exhibits, instead, a strong linear dependence in (pa)^, induced by the 0{a) contact term 

®The perturbative calculation of indicates explicitly that this coefficient is gauge-dependent. 
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Figure 2: Empty symbols denote S2l(p^), extracted directly from the lattice quark propa- 
gator. The filled symbols correspond to the subtracted form factor S2l(p^). Illustrated is 
the case with k, = 0.1344. 



proportional to c'g. This effect disappears in the subtracted form factor T,2l{p'^), defined 
from S~^{p), as shown in the same plot. 

As a further confirmation of the effectiveness of the subtraction, we show {l+z'^p'^)a2L{p'^)/p'^ 
in Fig. ^ for different values of n. At large and for all the values of the quark masses, this 
quantity has a good plateau corresponding, up to further C(a^) corrections, to the contact 
term ~ that we want to subtract. In the same figure, we also show o'2{p'^) / p^ ■, after the 
subtraction and extrapolated to the chiral limit (filled circles). This curve demonstrates the 
accuracy of the subtraction procedure, since in this limit we expect 'o^ip^^jp^ ~ {(lo) Ip'^- 



3 Numerical details and physical results 

In this section we briefly recall some elements of the lattice calculation which are explained 
in great detail in Ref. and present our physical results. These results are obtained 
on a sample of 200 quenched gauge fields configurations, on a 24'^ x 64 lattice and at 
(3 = 6.2. The value of the inverse lattice spacing, = 2.72(11) GeV, is obtained from the 
mx* mass. The quark propagators have been computed by using the non-perturbatively 
(!?(a)-improved Wilson action, for four different values of the light quark masses which 
correspond to the following set of the hopping parameters: 

K e {0.1352, 0.1349, 0.1344, 0.1333} . (27) 

All other details, concerning the analysis of the light hadron spectrum, can be found in 



Ref. [34 where a subset of 100 configuration was analyzed. 



9 



— 

1.5 - 



1.0 



0.0 



- • 
0.5 - • 



B ffiBB a B B B 5 B B B S B 5 5 



0.5 1.0 1.5 2.0 



(ap) 

Figure 3: The effect of subtraction of the contact term, c^, from the scalar part of the quark 
propagator. Squares denote (72/ p'^, as obtained from unsubtracted (albeit corrected by the 
(1 + z'^p'^) factor) propagators corresponding to the four values of the hopping parameters 
used in this study. The filled symbols denote d^2/p'^, obtained after the subtraction of the 
term due to 'S^ (see Eq. ( |26| ) ), and extrapolated to the chiral limit. 



As discussed in the previous section, Eq. (p6D involves the determination of the constant 
z. From the fit of our data to ( p^ in the interval 1.1 < {apY < 2, we obtain: 



z = {0.581(3), 0.591(2), 0.605(3), 0.639(2)} , (28) 

in decreasing order w.r.t. the K-parameter. With the value of z at hand, we follow the 
subtraction procedure described in the previous section and obtain the renormalized prop- 
agator, from which the quark masses can be derived. We now present the numerical results 
for Zq and for the quark masses in different schemes. 

3.1 The quark field renormalization constant Zq 

We fit the form factor Sii(p^) to a constant, in the large momentum region. According to 
Eq. (H), the value of that constant corresponds, up to tiny O^a"^) corrections |18|, to 



the value of the quark field renormalization constant, Zg, in the Rl-MOM scheme. From a 
fit in the region of large {apY G [1, 1.8] (corresponding to 7.4 ^ ^ 13.2 GeV in physical 
units), we get 

Zg = (0.849(2), 0.847(2), 0.846(3), 0.839(2)} . (29) 
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Figure 4: The lattice quark masses corresponding to the indicated k values, in the RI- 
MOM scheme. They have been obtained using Eq. (PT|). The error bars are not visible, 
because they are smaller than the symbols used in this Ggure. 



These values extrapolated (linearly) to the chiral limit give 0: 

Z(°)(2 GeV) = 0.853(3), (30) 

which is to be compared with zjf'^ = 0.875, from one-loop (boosted) perturbation the- 



ory Note that the quark mass dependence of Sii(p ), which according to Eq. 
comes from both the mass dependent term in Zg and the term proportional to cngi? is 
rather weak. 

3.2 Extracting the physical quark mass 

From the subtracted inverse propagator, S^^{p), we get the renormalized quark mass in 
the RI-MOM scheme which, according to Eq. (Rf), corresponds at large to 



m^^(/i) 



Si(p2 



(31) 



This quantity is unaffected by the correcting factor (1 + z^p^). In Fig. we display m^^(/i) 
for different values of k. The numerical value of m^^(/i) can be read off directly from Fig. ^. 
Note that m^^(/x) is derived in a completely non-perturbative way. 

^The value of Kcr^t is 0.135855(19). 
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In practice, to reduce the statistical fluctuations, we extract m^^(/io), from a fit in the 
interval 1.1 < (/ia)^ < 1.8, corresponding to 7.5 GeV^ ^ A*^ ^ 13GeV^, by using 



^^'(-") = ^WTT^ "^^'(/^o) , (32) 



with c^^(/i) computed in perturbation theory. We have chosen (/ioo)^ = 1.45, in the middle 
of the fitting interval, corresponding to /iq = 3.28 GeV. With this procedure the error 
induced by the use of perturbation theory is negligible, namely we find that the differences 
between NLO and N^LO are less than 1 %. The reason is that the expression in Eq. ( ^2]) 
depends on the ratio of c^^'s evaluated at different but close scales. 

The numerical values of c^^(/i) has been obtained by using the quenched expression for 
with Aqcd = 318 MeV i|. We checked that by using Aqcd = 238 MeV, as found 



in Ref. |T^, the central value of the masses is increased by less than 1%. 

To compute the physical values of the light and the strange quark masses, we invoke 
the procedure described in Ref. 0]. The masses are fitted as quadratic functions of the 
squared pseudoscalar meson masses. By using the method of "physical lattice planes" 0, 
we fix the average up-down and the strange quark masses, from the vr and K meson 
respectively The results are the following: 

(3.28 GeV) = 5.1(5) MeV , mf (3.28 GeV) = 118(9) MeV . (33) 

Since it is costumary to give the quark masses at the renormalization scale /i = 2 GeV, we 



use again Eq. (|32D to rescale m^^(3.28 GeV) to m^^(2 GeV). This time, we have to run the 
mass to a lower scale than the one used in the fitting procedure (7.5 GeV^ i$ A*^ ^ 13 GeV^). 
For this reason the uncertainty due to higher orders is larger than before, namely of the 
order of 4 %. The results are those given in Eq. (|^. We note, in passing, that there is 
no reason, if not for comparison with other calculations, to evolve the masses down to 
11 = 2 GeV. Indeed, it would be much more convenient to work at a scale /i > 3 GeV 
where the uncertainty induced by higher order perturbative corrections is negligible. 

We now illustrate the procedure adopted to obtain the quark masses in different schemes. 
The renormalization group invariant quark mass, which is a scheme and scale independent 
quantity, is related to m^^(/i) by the expression 



^ cRi(/x) ' ^ ^ 

whereas for the MS mass we have 

m^{^^) = c^^(/.) I = (/.) . (35) 



The resulting values of m^^^ should be flat in a large range of yU^, as confirmed by the data 
shown in Fig. ^. By using (^4]), we obtain the following results 

""^^Note that consistent results are obtained when the (^-meson is used to extract the strange quark mass. 
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Figure 5: Renormalization group invariant quark masses, obtained after dividing out the 
scale depending part in the RI-MOM scheme to N^LO accuracy, c^^(/i), from the quark 
masses depicted in Fig. ^. Dashed hnes correspond to the fitting interval used in ( |52D . 

NLO N^LO N^LO 

mf^^ = { 8.5(9); 7.8(8); 7.6(8)} MeV , 

mf^i = {196(15); 182(14); 177(14)} MeV . (36) 

In the MS case, from Eq. ( |33| ) we get the results quoted in Eq. (|TT]). 

In the calculation of m^*^^ and m^^(/i), we have used the physical value of ^^(/i), 
corresponding to as{Mz) = 0.118 computed with the appropriate number of active 
flavors, e.g. n/ = 4 at 2 GeV. This choice can be justified by assuming that the masses in 
Eq. (1^) are the physical ones, up to some unknown quenching errors. We checked, however, 
that the results by using the quenched with Aqcd = 318 MeV, would be different by 
less than 2 % in all the cases considered in this paper. 

Conclusion 

We have applied a new method to compute the renormalized quark masses from the lattice 
quark propagator using the OPE. We have discussed the subtleties related to the improve- 
ment of the propagator and especially the troubles arising from the presence of contact 
terms. Some of these problems could be avoided by working with Sl{x) instead than Sl{p)- 
Feasibility studies are underway. The main results, given in the introduction, are well com- 
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patible with values of the masses obtained by standard lattice methods 

^^(2 GeV) 



very good agreement with the recent result of Ref. ||3^ 



They are also in 
114 ±24 MeV, 

obtained by using the model independent QCD sum rule analysis at a N^LO (see also 
Ref. m). 
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Appendix 

In this appendix we list the formulae which have been used to compute the perturbative 
scale dependence of the quark masses. 

The effective QCD coupling is governed by the /^-function which is known to 4-loops: 



/i 



djj? 



TT 



n=0 



TT 



n+2 



(37) 



where the coefficients are 37 



* = J (ii -5"/ 1 . 



nMS 
P2 



oMS 
^3 



if 


2857 


64 V 


2 


1 


"149753 


256 


6 




/'50065 




I 162 



325 2 



3564C(3) - 



1078361 6508,, , , 
-^ + ^C(3))n,+ 



+ — ^C(3) ) n^f + ^n'f 
-1 ' ) f 729 ^ 



(3^ 



The coefficients of the mass anomalous dimension which describes the running of the quark 
mass, 



d 
d/j,' 



-m 



(„) f O^sifJ-) 



n=0 



TT 



n+1 



(39) 
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are also known up to four loops in both RI |T8[ and MS 
MS 1 /202 20 



schemes. We list them all: 



(Tm ) 
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(2)>MS 
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64 



1249 



3 g'^/j' M^)""' 
/2216 160,, A 



1 
li 
140 



52 

, 126 nf 

16 V 9 ^ 



-n. 



RI 



1 

64 



20911 3344,, , /18386 128,, A 928 ^ 



/91723 34192^,^, 18400^,^,, 

+ — C(3) - 880C(4) ^C(5) rif + 



V 27 
}h2A2 



9 



9 



^00 



160 



+ — C(3) - — C(4) ) - ( ^ - ^C(3) ) 



332 64 



(3)>M 



256 



300665987 15000871 , , 6160,, , 

<(3) + ^C(5) 



648 108 
/7535473 627127 



C(3) 



4160 



V 108 54 
/670948 _ 6416 



C(5) + 



3 

18832 o 
719 ^ 



243 27 

The corresponding evolution part in the running quark masses are: 



489 



242 



TT 



+ 



25335863 19 



1405536 6 



C(3) 



48247704573745 170324909 35' 
220410533376 2509056 ^ '* ^ 36^ 



TT 



(40) 
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TT 



(41) 



499 fa,{fi)\ 6375961 /a,(/i) 



726 



7r 



4216608 V vr 



344717507317 6293,,^, 25,,^, 

+ C(3) C 5 

55102633344 3564^^ ^ 6 ^ ^ 



for = 0, and 



TT 



+ 



5679460183 119 



337500000 



30 



TT 



C(3) 



(42) 



TC 



+ 



14533180260067051 1437607219 ^^/'(r:\ 
91125000000000 21600000 ' ^ T^*^ 



TT 



; (43) 
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/, N12/25 /i , ( O'sifi) 

1^ 3750 I TT 



+ 



793412683 
337500000 



:C(3) 



TT 



+ 



57222640693973 _ 2202791 ^ 7 ' 

7593750000000 ~ 337500 ^ 3^*^ ~ 18^^ 



TT 



} , (44) 



for Uf = 4. 
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